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■^ . The phase diagram of a simple area-preserving map, which was motivated by the 



quantum dynamics of cold atoms , is explored analytically and numerically. Periodic 



a, 



^ ' orbits of a given winding ratio are found to exist within wedge-shaped regions in 

f^ ■ the phase diagrams, which are analogous to the Arnol'd tongues which have been 

o, 

■^ , extensively studied for a variety of dynamical systems, mostly dissipative ones. A 

O 

^O I rich variety of bifurcations of various types are observed, as well as period doubling 

(~| ' cascades. Stability of periodic orbits is analyzed in detail. 

-(— > 

a • I. INTRODUCTION 

Periodic orbits are the backbone of the description of dynamical systems [l|,y,|3|, as they 
r\ \ enable global understanding of motion in the entire phase space. In the present paper we 

explore the periodic orbits of the map jj,|5|, |6|: 

Jn+i = Jn + k sin 9n+i + 27rf2 mod 27r 

6n+l = On + Jn mod 271 . (I. 1) 

This map will be denoted 9Jt in this paper. Like all the other maps which are introduced in 
this paper, it maps the 2-torus onto itself; however, "mod 27r" will always be left understood 
in all equations which follow. The map QJt reduces to the Standard Map J3[ for i7 = 0. When 
k = 0, for each rational value j/p of Q the map has periodic orbits of period p. For small 
k > such orbits still exist in a range of Q around j/p. This range of Q decreases with p. 
The region in the {k, Q) phase diagram, where periodic orbits of period p are found, has the 
shape of a wedge, with the tip on the Q axis {k = 0). This wedge grows wider as k increases, 



and intersects other wedges. Furthermore, as k increases more orbits are generated, others 
lose stabihty, and an extremely complicated {k, Q) phase diagram is generated (see Fig. 1 
of 1^). The purpose of the present paper is to understand this diagram. 
^ n the case of dissipative systems, this sort of phase diagrams have been studied extensively 
21, Isl, Is! • The standard sine-circle map p| : 

en+i = 9n - K sm{9n) - 27rQ, (1.2) 

is a paradigm in the study of mode-locking for dissipative maps, li K = and f2 is a rational 
number j/p any trajectory of the sine-circle map returns to its initial value (modulo 27r) 
after p iterations. For < i^ < 1, mode- locking is observed; that is, over a range of Q 
values around j/p (the mode-locking interval) a periodic trajectory with rational winding 
number j/p persists. This periodic orbit attracts all other orbits asymptotically in time, so 
that all of them eventually acquire this winding number. The widths of the mode-locking 
intervals are exponentially small in p, and increase with i^ up to i^ = 1. The regions thus 
formed in {K,Q) parameter space, terminating at K = 0, Q = j/p, are known as Arnol'd 
tongues [T], |9|. For K > 1 the tongues intersect, stability of orbits is lost, and new orbits 
are generated. In the same way as the sine-circle map p. 2|) is representative of dissipative 
systems, the map p. 1|) is representative of area preserving ones. The wedges in the phase 
space of p. Ij) where the periodic orbits of period p are found will again be referred to as 
Arnol'd tongues; however, it will be shown that the phase space structure, and the structure 
of the tongues in particular, is significantly different. This is expected, because orbits of 
area-preserving maps cannot attract. 

The map p. 1|) was obtained in studies of the quantum dynamics of cold atoms, which are 
kicked by a standing light wave, and are at the same time accelerated by gravity |5ll6llllLll2j|. 
The gravitational acceleration is proportional to Q, the strength of the kicking to k, and e 
is the separation from quantum resonance, which occurs when the ratio between the kicking 
period and a natural period of the center-of-mass motion of the atoms is an integer. It was 
shown that e plays the role of Planck's constant, and the limit e ^ (while k = ke is fixed) 
is a pseudo-classical limit that is useful to describe the dynamics for small e p, |6|, |lO| . In this 
limit the map (II. II) was obtained. Its periodic orbits, and the stable islands around them, 

nn 

account for the Quantum Accelerator Modes p, |a|, which were observed experimentally 

biyu- 



Compared to that of dissipative systems, the phase diagram of (II. 1|1 is much more com- 
phcated, mainly because overlaps of different tongues occur at arbitrarily small values of k. 
Therefore, its numerical exploration demands high precision in determining periodic orbits. 
A property of p. Ij) which proves of decisive help in this task is that p. Ijl . like the Standard 
Map, can be written as a product : 

M = 3b:^a , (I. 3) 

where the maps Jaj'^b are defined as: 

3a : Ue) -^ {-J,e + J) 

3b : UO) -^ {-J + ksm{e) + 27vn,e). (1.4) 

and are involutions, that is J^ = J^ = 1 (the identity map). If a map has property p. 3jl . 
then the search for its periodic orbits is greatly simplified, if restricted to the class of those 
orbits which include fixed points of either involution, as we review in detail in Appendix 1X1 
This method was applied in [l^ for the calculation of the last separating KAM torus of the 
Standard Map. In the present work this method is used to calculate periodic orbits of p. Ij) . 
both analytically, in the perturbative regime in Sec m and numerically in other regimes in 

secirm 



Many classic items of the theory of dynamical systems, such as period-doubling cascades, 
and other bifurcations |l|, |2, |3| are met in the study of map p. 1|) . In our analysis we have 
focused on tongues where Q at the tip (vertex) takes the values 1 (or 0), 1/2 and 1/3. 
We searched for periodic orbits such that one of their points is a fixed point of one of the 
involutions. As the boundaries of the tongue are crossed from inside, periodic orbits are 
annihilated as a result of coalescence of a stable orbit with an unstable one. Near the tip 
of the tongue j/p the orbits have j and p coprime. As k is increased these become unstable 
at some stability border where a period doubling bifurcation takes place. The ratio j/p 
does not change as a result of continuity. This period doubling is followed by a period 
doubling cascade, which satisfies a characteristic scaling rule of such a route to chaos for 
area preserving systems, and is universal [IS]- Also below the stability border a pitchfork 
bifurcation of an orbit of period p to two orbits of the same period takes place, resulting in a 
cusp in the {k, Vl) phase diagram (compare to [7|). The analysis is performed analytically for 
small k in SecOand mostly numerically in Sec. IIIII Examples of periodic orbits that cannot 



be calculated with the involution method are presented in Sec. IIVI The involution method 
is reviewed in Appendix |X] and various symmetries of the periodic orbits are presented in 
this Appendix as well as in Appendix El Properties of Gauss sums and the relation to 
points on periodic orbits are outlined in Appendix O The numerical method for finding 
periodic orbits is presented in Appendix O and a Hopf bifurcation is examined analytically 
in Appendix |E| As tongues intersect, there are regions in the {k,Q) phase diagram where 
orbits characterized by different j/p coexist. These orbits are found in different parts of the 
{9, J) phase space, and typically do not "interact". 

Here we are interested mainly in stable periodic orbits. If the issue of stability is ignored, 
strong rigorous results can be obtained. In particular, work by T. Oliker and B. Wajnryb 
jla] implies that if a line in the tongue (pi, ji) starts from its tip (vertex) and intersects the 
boundary of a tongue (^2? J2) which is its Farey neighbor [4,17-19], then it intersects along 
the way all the tongues {p,j) which are in the interval [—,—), intersecting all its Farey 
neighbors in the correct order. The same work also suggests that in the center of the tongue 
(j, p) all the orbits of period p ■ I {I integer) exist, and their winding number is j/p. 

II. PERTURBATIVE CALCULATION OF PERIODIC ORBITS 

In this section the periodic orbits of the map f ll. 1|) are calculated to the first order in k. 
Throughout this section J„ and 6n are understood to be on the 27r torus, notably they are 
taken mod 2tt] moreover, Q is also taken mod 1. First we note that the n— th iterate of the 
map ( II. Ijl at zeroth order in k is: 

Jn = Jo + ^TiQn 

n-l 

6I„ = 60 + ^ Jn' = 0o + Jon + vrfin(n - 1) . (II. 1) 

n'=0 

In order to determine J„ at the 1st order in k, only the 0— th order in 6 is required, and so, 
on substituting the 2nd eq. of ()II. 1|) in p. 1|) . and iterating n times, we find: 

n 

Jn = Jo + kY^ sin[^o + (^0 - 7rfi)n' + vrfin'^] + 27rfin. (II. 2) 

n'=l 

Substitution of this equation in the second equation of ( II. Hi yields 

n 

6)„ = 6*0 + Jon + 'nVLn{n - 1) + A; ^(n - n') sin[6'o + (Jo - T^^)n' + TrJln'^]. (II. 3) 



n'=l 



For a periodic orbit of period p the following equalities hold: 

Op -00 = 27TS 

Jp-Jo = 27TJ, {II. A) 

where the integers s and j are the winding numbers of 6 and J around the torus. At the 

0— th order in k, 

Q = ^- (II. 5) 

P 

must hold, and so periodic orbits can be found only for rational values of Q. For nonvanishing 
k, periodic orbits can be found in an interval around each rational value of Q. In what follows 
we solve for periodic orbits for small k, assuming that {Q — -) is also small (this assumption 
will be verified a posteriori) . Therefore we replace ^2 by - in the terms that are proportional 
to ~k in (ITOll and (HTHl . 

In order to find a periodic orbit of period p , we seek for a point {6q, Jq) which is at once 
a point in this orbit and a fixed point point of one of the involutions p. 4|) f Appendix IX|l . 
The fixed points of 3a satisfy 

Jo = (II. 6) 



while the fixed points 


oiJB 


satisfy 










Jo 


-¥ 


sin 6'o + 7rf2 


or 














Jo = 


-A; sin 


. 6^0 + vrf] + TT 



[II. 7) 



(II. 8) 

Since the periodic orbits of the map 971 of p. 1|) are identical to those of its inverse, we find 
in Appendix El that if in a periodic orbit of 371 there is a point with momentum Jj, then 
in a periodic orbit of this map, not necessarily the same one, there is also a point with 
momentum — Jj. If there is only one periodic orbit which is stable (or unstable), then the 
corollary at the end of Appendix [Bl implies that for each point with momentum Jj there is 
also one with momentum — Jj in this orbit, and so if p is odd then Jq = is in this periodic 
orbit. Therefore Jq of ()II. 6|) is in a periodic orbit of odd period p. On the other hand, for 
even p, ()II. 7|) or pi. 8|) is in the periodic orbit. This is discussed in Appendix 1X1 
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Froni( III. 2j) and ( III. 4j) we find that in the leading order in k (or e) a point on a periodic 
orbit of period p satisfies 



p lip 



e V 



[{Xp-l)n'+n'2] 



c.c. 



(II. 9) 



n'=l 



where Xp = fo^' odd p, where fill. 6|) holds, while Xp = 1 for even p where pi. 7|) holds and 
Xp = 1 + p/j where pi. 8|) holds. We will solve equation pi. 9|) for 6q to the first order in k 
and then show that it is consistent with pi. 3|) in this order. 

For this purpose it is instructive to use properties of Gauss sums 

consider 

p 



20 



21 



22 



23| . We 



(II. 10) 



which can be summed 



21 



,|2a|23, see Eq. (36) of j|) to give 



Q^ 



i(ij,p) 



G{p,j) = Vpe 



where ^{j,p) is real (see Appendix IH]) . Substitution in pi. 9jl leads to 



(II. 11) 



27r ( - - n 



Vp' 



■ sin "(9 



where 



^ = Oo + aj,p). 



(II. 12) 



(II. 13) 



A 6q that belongs to a periodic orbit must satisfy this equality. 

The Gauss sum pi. 10|) also satisfies (this is a corollary of Eqs. (34) and (35) of |^) 

p 



m=l 



^g«'^^[{Xp-l)m+"i2] 



^G{p,j) podd 



(II. 14) 



UG{p,j) + 1) P even 



Substitution in pi. 3jl for n = p, where p is odd, and using pi. 4|1 results in: 

k p — 1 



2lT 



p 2 



Vp 2 



■ sin ■(?. 



(II. 15) 



It is easy to see that it is consistent with pi. 12j) for s 



j(p-i) 



similar way that 



s n k^ . 

27r[---pJ = ^-sm'i? 



. For even p one finds in a 



(II. 16) 



(8,1)- 




(7,1) (6,1) (5,1) (4,1) P'2) (3,1)(8,3) (5.2) (7,3) (2.1) (7,4) (5,3) (*>'=> (3,2)(;,.j \^^ (5,4>(6,5) 



(8,7) 



FIG. 1: Various tongues, marked by {p,j), in the k, fi phase diagram. For small k the numerical 
results (dots) coincide with the analytically determined tongue boundaries (solid lines). 



and consistency with pi. 12|1 is possible for s = y. Equation ()11. 12|) determines the periodic 
orbits in small k. One finds solutions in the range 

~k 



\1 -Q\ < . 

P ~ 277^ ■ 

This region has the form of a wedge and for small k its boundaries are the lines 



(II. 17) 



Q 



J 



k 



± 
P 27r^ 



(II. 18) 



and they are presented in Fig. [T] 

Each tongue is characterized by p and j. The inequahty pi. 17|) justifies our assumption 
that I- — i7| is at most of order k (or e). We see that indeed for small k the boundaries are 



well described by fill. 18|1 . The point i) on the orbit is given by pi. 12|) . Special points are 
the two in the center of the tongue '(9 = 0, tt and the two on its boundaries '(9 = |, |7r. Since 
C,{j,p) is known, one can determine 6q from ^, as will be shown in Appendix O (see table 
In]). The specific values of ^o are of no importance for discussion in the rest of this section, 
since all the results may be expressed in terms of •&. In the regime of small k, for each value 
of k and fl there are two periodic orbits, one is stable and the other is unstable. The two 
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get closer and closer as the boundary is approached and coalesce at the tongue boundaries. 
We tested numerically the validity of ()II. 12|) for small k. 

The stability of an orbit can be determined as usual from the tangent map. The tangent 
map of p. Hi is 

(I + k cos 9 k cos 9 \ 
(11.19) 

and the tangent map over p iterates along a period— p orbit (also termed Monodromy matrix) 
is given by 

Tp{9o) = l[T{9,) . (11.20) 

A periodic orbit is stable if 

\Tr[Tp{9o)]\<2 (11.21) 

and unstable if the opposite inequality holds. In order to calculate the product pi. 20|1 to 
the first order in k we note that pi. 19jl can be decomposed as 

T{9,)=% + ~kTi{9,) (11.22) 

where 

- (^A 

r„ = (II. 23) 

and 

r,(,.) =(■=-"■ -^"'V (,,.24) 

With the help of the identities Tt{A + B) = TtA + TtB and Tt{AB) = Tt{BA) we find that, 
up to first order in k 

p 



TT{T,{9o)) = TTTl + ~kJ2'^T{Tr'T,{9,)). (11.25) 



-(P-1)7 

iri^j 
It is easy to see that 



p 1 
and 



n={^\ (II. 26) 



_ —(v-1) I p COS 9i cos 9i , 

T,{9,)T^^ '^=1'^ ' \. (II. 27) 





Using pi. 10|) . ()II. 11|) . and the fact that the angles along a periodic orbit are given to the 
zeroth order in k by the 2nd equation in ()II. 1|) . one finds that 

Tr(Tp(^o)) = 2 + kp^/^ cos^ (11. 28) 

where ■& is related to ^o by pi. 13|) . 

At the boundaries of the tongue given by pi. 18j) the stable and unstable orbits coalesce 
and Tr(Tp(6'o)) = 2. Inside the tongue the orbits with | < -i? < |7r are stable and the ones 
in the other part are unstable. This was numerically tested. 

III. NUMERICAL ANALYSIS OF PERIODIC ORBITS AND STABILITY 

BORDERS 

In the previous section the k, Q phase diagram was analyzed in the framework of pertur- 
bation theory for small k. This approach is valid near the tip of a tongue. In this section 
the phase diagram is studied in the non-perturbative region. In this region, the edges of 
a tongue deviate from the straight lines pi. ISjl (except in the case p = 1, where pi. 18jl 
is exact). Furthermore, this region is characterized by instabilities and bifurcations, which 
are observed as k is increased. We numerically explore this behavior for some low-period 
tongues. 

As periodic orbits coalesce at the edges of a tongue, their precise numerical determination 
is difficult near the edges; and this, in turn, makes the precise determination of the edges 
themselves somewhat delicate. An algorithm we have developed to locate the boundary is 
presented in Appendix O From ( III. 12j) one finds that if -t^s and "du belong to the stable 
and unstable orbits respectively, then, in the vicinity of the boundary, and for small k, 

i^u - ^s)' = levrVp^ (III. 1) 

k 

where 6Q is the distance from the boundary in Q. This equation can be used to check 
consistency of results. 

The region of existence of a stable {j/p) orbit is bounded on the left and on the right by 
the edges of the tongue, and from above by a stability border, which lies in the interior of 
the tongue ( see, e.g. Fig. El). Whereas at the edges of the (j/p) tongue, period p orbits, 
whether stable or unstable, disappear, at the stability border stable periodic orbits become 
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unstable, and, in the cases considered in this work, a period-doubhng cascade follows. 
Loss of stability is foreshadowed by the perturbative equation pi. 28|) . because the second 
term in it is negative for the stable orbit, and so it would lead to violation of pi. 2H) . and 
hence to instability, for \kp^^'^ cos{'d)\ > 4. According to this argument, the location of the 
stability boundary in k decreases with the period as p"^/^ (in agreement with earlier work 
^). However, this instability occurs for relatively large values of k, where perturbation 
theory is no longer valid, and therefore pi. 28|) may just provide hints there . In particular, 
it suggests that stable regions are more likely to be found close to the margins of a tongue, 
where cos('i9) is small. This is indeed observed in numerical results: as k increases, the 
stability border asymptotically approaches the edges of the tongue. 

All periodic orbits which are calculated in this section are found by the involution method, 
which is described in Appendix^ As explained in Appendix^ they are numerically found 
from the requirement that one point in the orbit is a fixed point of an involution. In the 
case of odd period p, this is involution 3a of p. 4j) : for even p, it is typically 3b of p. 4j) . 

A. The tongue j = 1, p = 1. 

For this tongue many results are known analytically, and were presented in p, la, |lO| . In 
particular the boundary of the tongue, and the stability border satisfy (p], Eq. (32)): 

h = 27r|fi-l| , (III. 2) 



k, = V16 + 47r2(fi-l)2 . (III. 3) 

These are shown in Fig. |21 along with numerically obtained data (the latter are shown 
in order to check consistency of our numerical method of finding the boundary). As k 
increases beyond kg a period doubling bifurcation takes place, such that the period 1 orbit 
turns unstable and a period 2 stable periodic orbit appears. Around the p = 2 orbit a 
resonant orbit of period 2 ■ 3 is formed, as is shown in Fig. |3] for the case marked by a 
circle in FigEl This mechanism is discussed in Appendix |El For Q = 1, k = 27t a. pitchfork 
bifurcation of the p = 2 orbit into two p = 2 orbits takes place. The regions where such 
orbits exist are bound by dashed lines in FigO They are magnified in Fig. El and described 
in detail in the corresponding caption. 

The bifurcation which is observed on crossing the stability border upwards is actually the 
1st in a period-doubling cascade (marked in Fig. |21at i7 = 1.3692 by rhombi). Let us denote 
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by 56n the difference in 6 between the points separated by half a period in the orbit of period 
2"p, when it becomes unstable; p is the original period at the start of the period-doubling 
cascade. We follow the cascade up to n = 4 and find that the numerical results as n increases 
are consistent with: 

~kn-~koo-^ (HI. 4) 

sen ~ \ (III. 5) 

where kn are values of k where the {n + 1)— th bifurcation is observed. Estimates for 5 
and a are reported in Table HJ which summarizes this sort of estimates for other tongues 
as well, fcoo was calculated using the ratio between the LHS of pil. 4|) for n with the one 
for n + 1, thus eliminating the proportionality coefficient. The scaling factors 5 and a are 
universal numbers for the period-doubling route to chaos, in area preserving systems with 2- 
dimensional phase space. The theoretical prediction of the renormalization group treatment 
is 5 ^8.721 a^ -4.018 [15|. 

B. The tongue j = 1 ,p = 2. 

The structure of this tongue is presented in Figs. El and IHl First we study the phase 
portrait at the center of the tongue {VL = 1/2) for various values of k. For k = 3.0 (marked 
by "o" in Fig. EJ, which is below the stability line , the phase portrait is presented in Fig. 
m (2nd and 3d slice from top). In the upper slice, the function Fe^Oo) is plotted, the zeros 
of which are used to locate 6o of periodic orbits according to the Involution method, as 
explained in Appendix O A similar analysis was performed for all periodic orbits. As we 
increase k, but keep it lower than the stability border ks, a new stable orbit of period two 
appears at A; > vr. Further bifurcations are found for k > tc, and these are described in Fig. 

El 

The bifurcation leading from one stable period-2 orbit to two distinct stable orbits with 
period 2, which is numerically observed at fc = vr, is easily obtained analytically. It is a 
standard pitchfork bifurcation. To see this, we iterate the map ()I. 1|) twice. The condition 
for a period two orbit (modulo 2n) is 

k sin 6*0 + A; sin ( 610 + - sin 6*0 + - = 0, (III. 6) 
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p=i i=i 




FIG. 2: Phase diagram of the (j = l,p = 1) tongue. Full lines show the analytical edges of 
the tongue (straight lines emanating from the point marked by 1 on the horizontal axis) and the 
analytical stability border of period 1 orbits. Superimposed dots show numerical results. The 
meaning of symbols (circle, rhombi) is explained in the text. The dash-dotted line marks the 
stability border of the p = 2 orbit which emerges at the period doubling of the p = 1 orbit. The 
region delimited by the dotted-line rectangle is described in Fig. El 



where fl = 1/2 and pi. 7|) were used. When k < it there is one stable trajectory. At 
k = 71, this trajectory is {(7r/2,7r) , (37r/2,7r)} and is marginally stable (the trace of the 
linearized map Tr[T2(6'o)] = 2). When k = it + 6k, two new stable trajectories emerge. To 
find these zeros of Eq. pil. 6j) we replace 6o = vr/2 + 66o and neglect terms of order of 



0{6k69o 



0{59q'^) and higher. This yields the following cubic equation for SOq: 



S9o 




0. 



(III. 7) 



The same equation is obtained from the substitution Oq = 37r/2 + SOq. Eq. pil. 7|) yields 
three periodic orbits with period 2. One of them has 69o = 0. It continues the unique 
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FIG. 3: A zoom-in of the region delimited by the rectangle in Fig. |2 Stable period 2 orbits are 
found below the line JHG and inside the loop BHB, as well as between the lines EF and BF', and 
between BD' and CD. Inside the kite-shaped region AEBC, two stable orbits of period 2 are found. 
A very complicated behaviour is found in the vicinity of H where period 2 and period 4 stable 
orbits coexist. 



orbit at k = TT and, like that orbit, it is marginally stable at this order, numerically, it was 
found that it is actually unstable, and that ^ 59q = o{6k). The other two orbits have 
S9q = ±W ^. They are stable at this order, sharing the same value of Tr[T2] = 2 — 27r5k, 
and are numerically found to be really stable. 

On further moving upwards along the axis of the tongue ( that is, increasing k while Q = 
1/2), the two period-2 stable orbits persist until point B is reached. There they turn 
unstable, and give rise to two stable period-4 orbits. Two stable period-2 orbits (which we 
term "a" and "6" for convenience) are actually observed inside the whole kite-shaped region 
AEBC in FiglHl Each of them is observed , alone, in a larger region. For instance , "a" 
is observed in between the lines AD, ED' and "6" is is observed inbetween AE and CE' . 
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FIG. 4: Phase portrait for k = 5.5825 > kg, and $7 = 1. 



On the hne AD , "a" coalesces with an unstable orbit and disappears, and "6" disappears 
in the same way along the line AF. Across the line ED', "a" bifurcates to a period-4 orbit, 
and the same happens to "6" across CF' . Resonances of order 3 similar to those found for 
the p = 1 tongue and are presented in Fig. |3] are found also here . As fc is increased above 
the stability border (dot-dashed line in Fig. [SJ, a period doubling cascade is started. This 
was studied in the same way as for the p = 1 tongue. Results are shown in Table U and are 
marked by rhombi and triangles in Figs. ElandlHl 

So far the periodic orbits in the tongue p = 2, j = 1 were found identifying one point of 
the orbit as a fixed point of 3b- There is also at least one stable periodic orbit that none 
of its points are fixed points of 3b, but some are fixed points of 3a- For i7 = 1/2 one such 
an orbit of period 4 is {Xq = (0, 0) , Xi = (vr, 0) , X2 = (0, vr) , X3 = (vr, n)}, as is easily 
checked . The points Xq and X2 are fixed points of 3a while Xi and X3 are interchanged 
by this involution. All this demonstrates the Corollary in Appendix El We checked (by 
calculating the trace of the tangent map) that the orbit is stable for < A; < 0.7320 and 
2.7320 < A; < 2.8284. 
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FIG. 5: Phase diagram for the p = 2, j = 1 tongue. The region of coexistence of two stable 
periodic orbits of period 2 is framed by the dotted-hne rectangle, and zoomed in Fig. IHl The 
thick dashed lines mark the tongue boundaries (found numerically), and the dash-dotted line is 
the stability border (found numerically). The solid lines indicate the perturbative analytical result 
of Eq. (|11. 18(1 for the boundary of the tongue. The small rhombi and triangles indicate the period 
doubling cascade for 0, = 0.6376 and 0, = 0.3624 correspondingly. 

C. The tongue p = 3, j = 1 



Our numerical findings about the p = 3, j = 1 tongue are summarized by Figs. IHlandlHl 
The stable period-3 orbit is lost going upwards across the stability border, the details 
of which are shown in Figs. |H1 and IHl . However, a stable such orbit reappears on further 
increasing k; it is observed in the region bound by the lines LR and L'R'. Between A and B 
in Fig. El there exist two stable periodic orbits of period 3 which we call "a" and "6". Both 
turn unstable at B. "a" coalesces with an unstable orbit and disappears on the line AD, 
and "6" disappears in the same way along AF. When crossing CF', "6" loses stability, and 
the same happens to "a" along ED'. 
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FIG. 6: A zoom-in of the region delimited by the dashed-hne rectangle of FigEI At the dash- 
dotted lines period-2 orbits lose stability, and stable period-4 orbits are formed. On the dashed 
lines coalescence of one of the period-2 orbits with an unstable periodic orbit occurs. This stability 
loss at the dash-dotted lines is followed by a period doubling cascade. The small rhombi and 
triangles indicate this cascade for i7 = 0.5122 and fi = 0.4878 correspondingly. 

IV. ORBITS THAT ARE NOT RELATED TO INVOLUTIONS. 



In the previous sections periodic orbits that include fixed points of 3a or 3b were found. 
An obvious question is whether periodic orbits can be found, that do not include such fixed 
points at all. First we find such orbits for fc = 0. For fi = j'/p, a point ( Jq, ^o) lies in a 
periodic orbit of period p if Jop + vrj(p — 1) = 27rs with s integer, and ^o is arbitrary, as one 
can see from p. 1|) and pi. 3|) . It is easy to find numerically values oi j,p, s such that none 
of the points in such an orbit is fixed under 3a or 3b'- for instance, with p = 6, j = 3, and 
s = 2, it is easy to check that none of the fixed-point conditions pi. 6|) . pi. 7|) or pi. 8|) is 
satisfied for any of the points in the periodic orbit. Such an orbit can then be used to find 
periodic orbits where none of the points is a fixed point of the involutions, also for A; > 0. 
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FIG. 7: Phase portrait for p = 2, j = 1, k = 3.0 and Q = 1/2. The upper shce presents the 
function Fq{9q) . The middle shce shows the period two stable orbit and the lower slice is a zoom 
on one of its island. 



For this purpose take A; > but small. Iterate p times the line of the fixed points of the 
map with fc = 0, generating a new line M^( Jq, ^o)- Let ( Jq, ^o) be a point of the intersection 
of these two lines. Although this is not necessarily a fixed point of M^ for /c > 0, it is 
reasonable to look for one in its vicinity, and indeed we found such orbits of period p = Q. 
They are plotted in Fig. [HI These orbits are also related by the symmetry J — > —J 
discussed in Appendix |B| 
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FIG. 8: Phase diagram for the p = 3, j = 1 tongue, indicating the region of coexistence of two 
stable period-3 orbits, which is zoomed in Fig. |^ The soUd line marks the tongue boundaries found 
numerically, and the dash-dotted line is the stability border also found numerically. The dashed 
lines bound the region in which two period-3 orbits coexist. As each of these lines is crossed from 
the inside of this region, one periodic orbit coalesces with an unstable orbit, and is annihilated. 
At the stability border a period doubling cascade begins. The rhombi and triangles illustrate such 
a cascade for il = 0.3936 and il = 0.2285 correspondingly. Between lines LR and L' H a stable 
period 3 orbit is again observed. 

V. SUMMARY AND DISCUSSION 



In this paper the tongues with jiv = 1) 1/2, 1/3 were analyzed in both the perturbative 
and the nonperturbative regime. A large variety of scenarios for the disappearance of orbits, 
the loss of their stability and bifurcations was found. We expect similar behavior for the 
higher order tongues, but the verification is left for further research. For higher values of k 
some novel results may be found but such an exploration is out of the scope of the present 
paper. 
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FIG. 9: A zoom-in of the region delimited by the rectangle in Fig. |HJ 



tongue (p,j) 


^ 


kg 


"^oo 


a 


5 


IT'max 


(1,0) 


0.3692 


4.6239 


5.1290 


-4.0645 


8.5839 


4 


(2,1) 


0.3624 


2.0019 


2.2790 


-3.9919 


8.8307 


4 


(2,1) 


0.6376 


2.0019 


2.2790 


-3.9919 


8.8307 


4 


(2,1) 


0.4878 


3.9126 


3.9711 


-4.0339 


8.7034 


4 


(2,1) 


0.5122 


3.9126 


3.9711 


-4.0339 


8.7034 


4 


(3,1) 


0.3936 


1.03 


1.1601 


-4.1375 


8.8534 


4 


(3,1) 


0.2258 


1.5169 


1.6167 


-4.0315 


9.0737 


4 



TABLE I: As mentioned in the text, the stability of orbits is lost through a period doubling cascade. 
In this table some examples are given ,for different values of 0,. 



The research in the present work was confined to the classical regime. The fingerprints 
of the rich classical behavior that was found on the quantum dynamics in the semiclassical 
regime may be of great interest. In particular the effect of bifurcations on the eigenstates 
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Q=0.5 k=0.345 




FIG. 10: The two period-6 orbits found at fi = 1/2 for k = 0.345. The two hnes indicate the 
invariant hnes under the Jb involution (as given by ()II. 7|) and (|II. 8|) ). thus exhibiting that none 
of the points of these periodic orbits is a fixed point of the 3b involution. It is also evident that 
none of the points is a fixed point of 3a, since none of the points is located on the invariant line 
of 3a which is given by (|11. 6|) . Points in one periodic orbit are marked by aj i = 1,2, ...,6 and 
points in the other by 6j. The involution 3a transforms the former orbit into the latter, as follows: 
ai — > bi, 02 — > ba, 0,3 — > 65, 04 — > 64, 05 — > 63, a^ — > b2, while 3b transforms them as 
follows: ai — > 62, 0,2 — > bi, 03 — > be, 04 — > 65, 05 — > 64, uq — > 63. This is a manifestation of 
theorem El of Appendix 



or the scattering states of the quantum system may reveal novel results. Effects similar to 



the ones found for quantum ratchets |2^ and "flooding" [25'] are expected. Also states on 
islands around periodic orbits related by symmetry may result in states of novel properties. 
These may lead to effects observable in atom optics experiments. 
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APPENDIX A: THE INVOLUTION METHOD 

This method allows for establishing an analytical relation between the coordinates J, 9 of 
a point on a periodic orbit, and in this way it considerably simplifies the task of computing 
periodic orbits. The method was introduced by Greene [1J| who implemented it in the 
computation of the approximants of the last bounding KAM trajectory of the Standard Map. 
In this Appendix the method is reviewed for the reader's convenience, and its application 
to the map OJt of p. Ij) is formulated. 

Consider a map that can be written as a product of two involutions, i.e. 

971 = ^2 2fi , (A. 1) 

where the maps Ji and 3^2 satisfy 31 = 3^ = 1 (the identity map). The following theorems 
hold: 

Theorem 1 //, for some point Xq in phase space, and some integer N, 

3, (Xo) = Xo , (A. 2) 

and 

3, (mt^ (Xo)) = on^ (Xo) , (A. 3) 

then 

mt^^ (Xo) = Xo . (A. 4) 
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Theorem 2 //, for some point Xq in phase space, and some integer N, 

32 (Xo) = Xo , (A. 5) 

and 

3, (971^ (Xo)) = OJl^ (Xo) , (A. 6) 

then 

gjt^^+i (Xo) = Xo . (A. 7) 

These theorems are easily proven by exphcit use of (jA. 2|1 and ()A. 3j) in the RHS of ()A. 4jl . 
and of dHH), and of (IX~6l) in the RHS of (IXTTj) . 

The opposite direction is not correct. Orbits can be found, such that no point in them is 
fixed under 3i or 32 (see section HVT) . 



Theorem 3 Let O = {Xi,X2, ....Xp} be a primitive periodic orbit of period p of the map 
971. Fori = 1,2 denote O^"^^ = 3i{0). Then O^"^^ is a primitive periodic orbit of period p, 
andO^^^ = 0(2). 

Proof From ^K~l^ it foUows that 971^2 = ^2971"^ Therefore, for any Xj E O, 

97lp2(X,)) = 32{m-\X,)) e 32{0) = 0(2) , (A. 8) 

hence 9Jl(0(2)) c 0^'^\ and then 9Jt(C(2)) = (9(2) follows, because 9Jt is one-to-one. Further- 
more, ()A. 8|) entails that O'-^^ is the orbit of any of the points in it, because this is true of 
O. Finally, 

OW = 3i{0) = 3i3232{0) 

= 3i32{0^^^) = Tl-\0^^^) 

= C(2) . (A. 9) 

Theorem IHl has the following direct 

Corollary: If in O there is a fixed point of either involution, then O = O^^^ = 0^'^\ and so 
the orbit is invariant under both involutions. The orbit can be divided in classes such that 
each class either contains one fixed point of one involution or it contains two points which 
are mapped onto each other by that involution. In the case when the involution is 3a, points 
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in each such pair have opposite J . For odd p, this imphes that a point with J = belongs 
to the periodic orbit (and is a fixed point of 3a) ■ 

Theorem ^ is employed for finding periodic orbits of even period and Theorem El for 



finding periodic orbits of odd period [l^. One solves for periodic orbits, assuming that in 
a periodic orbit there is a fixed point Xq of either 3 a or J^ (in our case, Xq = (Jo,6'o)). 
This assumption reduces the number of unknowns. Since the periodic orbits of QJt and of 
0K~^ are identical, and 2Jt~^ = 2J1J2 holds, we can interchange 3a and 3b in the calculation 
of periodic orbits. For orbits of odd period one can easily check that Xq = ( Jq = 0, 6*0) 
(eq. III. 6jl is a fixed point of 3a- In our applications periodic orbits with an odd period p 
contain a fixed point of 3a where pi. 6|1 is satisfied, while periodic orbits with an even period 
typically contain a fixed point of 3b, where (III. 7|1 or pi. 8j) is satisfied. By the corollary, 
these orbits must contain more than just 1 fixed point of the latter involution. 

APPENDIX B: TIME-REVERSAL SYMMETRY. 

In this appendix a symmetry of periodic orbits of OK is discussed, which is related to 
time-reversal. The map DJl may be written as 9Jt = DJIk'^f, where 

Mf : (J, 9) -^ {J, 9 + J) , 

TIk : UO) -^ {J + ksm{e) + 27Tn,e) . (B. 1) 

(throughout this appendix, mod 27r is left understood in maps). Let OJl' = DJIf^k- The 
map DJl' is easily seen to satisfy: 

m~^ = 3m' 3 

where 

3 : {J, 9) ^ {-J, 9). (B. 2) 

and so, denoting !H = 3DJIf, the map fH is an involution: IH^ = 1, and satisfies 

m-^ = D^mm. (B. 3) 

Therefore DK achieves the time-reversal of the map DJl. We have thus proven that: 

Proposition 1 971 zs time-reversal invariant, in the sense that an area-preserving involution 
91 exists, such that 9^9711H = VJl^^. Therefore, if O is a periodic orbit of VJl, then the 
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time-reversed orbit 9l((9) is also a periodic orbit of DJl, and the two orbits have the same 
stability. 

Corollary: If for the map p. Hi there is only one periodic orbit of a given stability (as 
determined by the tangent map), then O = 91(0) must hold, and so for each point in the 
orbit with momentum J there is another point with the value —J. 

APPENDIX C: GAUSS SUMS AND PERIODIC ORBITS 



In this appendix some properties of Gauss sums that are relevant for the present paper 
will be reviewed. Specifically we focus on the sum pi. 10|1 . that is related to 

p 



GiJ,P^l) = lT.^ 



- [Im+jm 



(C. 1) 



m=l 



where in our case / = j{Xp ~~ 1) vanishes for even p if pi. 7|) holds and it is equal to p if 
pi. 8|) holds while it takes the value —j for odd p. Using elementary calculations one can 
show 1^ that in our case ()C 1|) is just pi. 11|) divided by p, namely 



G{J,P) =pGU,p,l). 



(C.2) 



The calculation of ^{j,p) is more complicated and will be summarized in what follows. In 

I), sums like fid 1|1 are calculated for 



the standard number theoretical literature 



I = O^eve. , aud odd p^ Haunay aud Berry M gerreralized to other cases of Q (See 

also ^21] ). They show that this sum reduces to the average of the summands and is given, 
for j and p mutually prime, by Eq. (14) there, which reads: 



1 



G{j,p,l)= < 



Vp 



Vp 



p 



p 



p 



exp 



-fb-i; 



2/ l\2 



exp 



'-fij\p) V, 



p odd, j even, / even 



-f (p - 1) - '-^{4j\p)H^ p odd, J odd, / odd (C. 3) 



exp 



f J - 'fUWii? 



p even, j odd, / even 
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where {j\p) is the inverse (mod p) of j, notably 

j{j\p) = 1 mod p (C. 4) 

which can be calculated with the help of (whenever j is mutually prime with p) : 

{j\p) = f(p)-^ mod p (C. 5) 

where <^(p) is Euler's totient function, that is, the number of integers less than p that have 

(A 

no common devisors with p. \ is the Jacobi symbol, where a and b are integers so 

vJ 

that b is odd and does not divide a. The Jacobi symbol assumes the values ±1 only, and is 

(A 

the product of the Legendre symbols I for all the prime factors q of b. The Legendre 

Vv 

symbol is defined by 

a \ 1 if there is an integer m such that m^ = a mod a 

\^{ '^ (C. 6) 

q I I —1 if there is no integer m such that m = a mod q 

A more straightforward way of determining the value of a Legendre symbol is given by the 
next formula, which is based on Fermat's little theorem : 

= j^~ mod p (C. 7) 

P J 

keeping in mind that p is an odd prime and that j is not divisible by p. For our purposes it 
is convenient to write the Jacobi symbol in the form 

(C. 8) 





It is useful to note that 



and that 



and also that 




(C.9) 



(C. 10) 



(C. 11) 
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for all integers j. Other useful formulae are given in appendix B of J20| and in 22 1 and 23| 
(see also |2l|). 

We note that in our case, for even p, / = in case pi. 7jl . and / = p in case pi. 8j) . while 
for odd p, I = —j. In this case HC Hi takes the form 

1 



G{j,p,l) = — e^^^^'^^ 



(C. 12) 



where, for even p corresponding to pi. 7|) . 



ap,j] 



and for even p corresponding to pi. 8|) 



TT 



. 2 



(C. 13) 



^(P>j) = 77 



V 



p 



(C. 14) 



while for odd p 



ap,j] 



J 
p 



|(P-1)-|J(J\P)^^^ 



Up - 1) ~ 'ii^jWi' 



J even 



(C. 15) 



P 



j odd 



The phase ^{p,j) determines the relation between 9o and {} via pi. 13j) and can be used to 
check the accuracy of the calculations. At the center of the tongue ■(? = or -(9 = vr while at 
the boundaries {} = ^ or {} = ^. In table |n] we compare the values of the numerical results 
A6 = d — 9q for k = 0.001 and 6o at the center of some tongues with ^{p,j) calculated from 

(irrT3|) and dnnsi). 

APPENDIX D: A NUMERICAL METHOD FOR THE CALCULATION OF 

PERIODIC ORBITS 



In the numerical solution for periodic orbits we make use of pi. 6j) . pi. 7j) and pi. 8|1 in 
order to eliminate Jo- Then we impose the periodic orbit conditions pi. 4j) and require that 
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tongue {p,j) 


A9 = ^-eo 


e(p,i) 


(2,1) 


0.7855 


tt/A w 0.7853 


(3,1) 


0.5235 


it/6 « 0.5235 


(3,2) 


-0.5235 


-^/6 ^ -0.5235 


(4,1) 


0.7855 


7r/4 w 0.7853 


(4,3) 


2.3567 


37r/4 « 2.3561 


(5,1) 


0.6284 


7r/5 « 0.6283 


(5,2) 


1.2573 


27r/5 « 1.2566 


(5,3) 


-1.2573 


-27r/5 ^ -1.2566 


(5,4) 


-0.6284 


-tt/5 ^ -0.6283 


(6,1) 


0.7857 


7r/4 « 0.7853 


(6,5) 


-2.3558 


-3^/4 PS -2.3561 


(7,1) 


0.6733 


37r/14 « 0.6731 


(7,2) 


-0.2251 


-7r/14 w -0.2243 


(7,3) 


2.0209 


97r/14 Ri 2.0195 


(7,4) 


-2.0209 


-97r/14 PS -2.0195 


(7,5) 


0.2251 


7r/14 « 0.2243 


(7,6) 


-0.6733 


-37r/14 PS -0.6731 



TABLE II: Numerically computed values of A^ = 'd — 0q, compared to CiP^J) (see (|II. 13|) ). 



6'o is a common zero of the functions 

p 
Fj(^o) = A; ^ sin ^„ + 2'KQp - 27rj, 

n=l 

where j is determined by the tongue f ()II. 4|) and ()II. 5|) ). and 

p 
-Pe(^o) = pJo + k 2J(P ~ '^) sin 6'„ + 7iQp{p - 1) - 27rs 



(D. 1) 



(D.2) 



n=l 



that are obtained after iterating the map p. 1|) p times. We solved for zeros of Fj{6q) 
and sin[|F6)(^o)]- This is because j is known, while s is unknown and since the function 
Fgi^Oo) should be taken mod 27r. Then only the common zeros of Fj{9q) and sin[|F6)(^o)] are 
kept. The elimination of Jq by the involution method, described in Appendix ^ reduced 
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the problem to one of finding zeroes of a function of a single varia ble, that is a relatively 
easy task. We used the MATLAB7 routine "fzero" for this purpose [26]. In order to locate 
the boundary, we note that the zeroes corresponding to the stable and unstable orbits 
are separated by an extrenium of the function in question, say Fj{9q). At the boundary 
these zeroes coalesce and at this extremum the relevant function vanishes (at the extremum 
Fj = 0). The extremum can be easily found by the MATLAB7 routine "fminbnd" J26!]. The 
boundary of a tongue in the fc, VL space is the locus of points of vanishing extremum. It was 
found by successive change of k and Vt starting from vertex A; = 0, Q = j/p. The boundary 
was found by the overshoot method, making use of the fact that the value of the function 
at the extremum changes sign when the boundary of a tongue is crossed. 

APPENDIX E: THE RESONANCE OF ORDER 3 FOR THE p = 1 TONGUE 

When the parameter k is increased across the stability border, the typical scenario of a 
Hopf bifurcation is observed, with period doubling accompanied by the birth of a resonance 
of order 3, as shown in Fig. |3J Here we obtain analytical expressions for the first period 
doubling, and the normal form for the resonance of order 3. 

For the map p. Ij) that is studied in the present work, the fixed point ( Jq, ^o), discussed 
in subsection IIII A| loses stability for k^ cos 6*0 = —4. In order to analyze the resulting 
bifurcation, we study the small deviations 

69 = 9-00 

5J=J-Jo. (E. 1) 

From p. Ij) we find that to the second order in 69, 

6Jn+l = 5Jn + A69n + D{69nf 

S9n+i = 69^ + 6Jn+i , (E. 2) 

where A = kcos9o, D = —(k/2)sin9o. It is convenient to introduce the variable x = 

^ , n 

{D/2)69 and turn Eqs. ()E. 2jl into one second order standard equation (see j2|, p. 243) 

Xn+i + Xn-i = 2Cx„ + 2x„ , (E. 3) 

where C = {2 + A)/2. In our case we are interested in the vicinity oi A = —4, or C = —1. 
For C > —1, the fixed point Xn = (corresponding to the original point {Jo,9o)) is stable. 
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For C < — 1, it becomes unstable and a stable periodic orbit which consists of the points 

X2± = a±b, (E. 4) 

where a = -(C + l)/2 and b = |a/(C+ 1)(C-3) appears. If C = -I - 1], then a = r]/2 
and b = ^ \/'f]^ + ^)- To study the small deviations, we denote 

Xn = X2± + AXn , (E. 5) 

where we numerate the steps so that for even n the point Xn is near X2+, while for odd n it 
is near X2-. Iterating the map and keeping terms to the second order in the Axn, one finds 
for even n 

qn+2 + qn-2 = 2C' Qn + 2ql , (E. 6) 

where 

qn = 2[{C + 2X2-) + {C + X2^f]AXn (E. 7) 

and C = —2C^ + AC + 7. In our case C" = 1 — Sr/ — 27f. For odd n the same result is found, 
but with X2~ and X2+ interchanged. 

The map ()E. 6|) can be obtained from the Hamiltonian 

^2 



H = ^- + V{q) J2 ^(t-^) (E-8) 

n=— 00 

with the time rescaled by the factor of 2 and 

V{q) = -^A'q' + ^q' (E. 9) 

with A' related to C via C = {A' + 2)/2. To see this, one notes that the map obtained from 
the Hamiltonian ()E. 8j) is similar to Eq. ()E. 2j) . and therefore, Eq. ()E. 6|1 can be obtained 
from Eq. ()E. 8j) in the same way as Eq. ()E. 3j) is obtained from Eq. ()E. 2|) . The parameter 
that controls the potential is 

A' = 2C' -2 = -4C^ + 8C + 5 = -16r] - Arf . (E. 10) 

Using the Poisson summation formula we write the periodic 5 function in the form 
Xl^-oo ^it — n) = 1 + 2 X]^i cos(27mt). For small q the averaged, or time independent. 



part reduces to a harmonic oscillator with the frequency u = \/—A' = 2\/^r] + rf. The 
third order resonance is found for ui = 27r/3. Introducing the action-angle variables (/, 0) 
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for the linear oscillator part, where q = -^21 juj cos <\) and p = \/2Iuj sin </>, the Hamiltonian 

fp~8|l reads 

4 / J \ 3/2 
i7 = a;/+-(— j cos(30 - 27rt) + 7^(/, 0, t) , (E. 11) 

where 7^(/, 0, t) is the non-resonant part of the Hamiltonian ()E. 8|) . Introducing the resonant 
phase ih = (h — 2TTt/3 and the detuniner from the linear resonance S = a; — 27r/3, and denoting 
r = (4/3)/(2c<j)3/2 we obtain the normal form for the 3-resonance |27Ll28l|: 

H3 = 6-I + rl^'"^ cos(3?/') . (E. 12) 

Denoting x = v/cos(^), y = v/sin(?/^), we obtain another well known expression for the 
normal form of the resonance of the order of 3 (2 71 . l28l | 

H^ = S{x'^ + f) + r{x^ - 3xf) . (E. 13) 



As k increases the orbit (jE. 4|) loses its stability and bifurcates, typically via a period dou- 
bling bifurcation cascade. 
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